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Abstract 

For the Ising model, the calculation of the spontaneous magnetization leads to the 
problem of evaluating a determinant. Yang did this by calculating the eigenvalues in the 
large-lattice limit. Montroll, Potts and Ward expressed it as a Toeplitz determinant 
and used Szego's theorem: this is almost certainly the route originally travelled by 
Onsager. For the corresponding problem in the superintegrable chiral Potts model, 
neither approach appears to work: here we show that the determinant Dpq can be 
expressed as that of a product of two Cauchy-like matrices. One can then use the 
elementary exact formula for the Cauchy determinant. One of course regains the known 
result, originally conjectured in 1989. 
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1 Introduction 

Extrapolating from the Ising case and from series expansions, Albertini et al conjectured 
in 1989 Q] that the order parameter or spontaneous magnetization of the solvable N- 
state chiral Potts model is 

Mr = (1 - fc'Y (JV - r)/2iV2 - (1.1) 

Here r is an integer, between and N, and k' is a parameter that is "universal" 
in the sense that it is the same for all rows and columns of the lattice, even for an 
inhomogeneous model where the rapidities vary from row to row and from column to 
column. [2] Here we consider the ferromagnetically ordered phase of the system, where 
< k' < 1. This k! is small at low temperatures (high order), and tends to one at the 
critical temperature (vanishing order), so we can regard it as a "temperature". 

The author was able to derive the result (|1.1[) in 2005, 3 , 4 using an analytic method 
based on functional relations satisfied by generalized order parameters in the large- 
lattice limit. 
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If the vertical rapidities of the homogeneous model take a particular value, or if the 
those of a model with alternating vertical rapidities satisfy a certain relation, 5 , 6 then 
we obtain the "superintegrable" case of the chiral Potts model. For this case, Gehlen 
and Rittenberg 7 showed that the horizontal and vertical components of the transfer 
matrix satisfy the "Dolan-Grady" condition, 8 which ensures that they generate the 
Onsager algebra. This is the algebra generated by the transfer matrices of the Ising 
model. eqs. 59-61], [101 eq. 4.12], [11] - [18] Onsager used it to calculate the free 
energy of the Ising model. 

In fact the superintegrable chiral Potts model looks very much like the Ising model. 
(For N = 2 it is the Ising model.) If one considers the model on a cylinder of L columns, 
with the spins on the top and bottom rows fixed to the value zero, then the row-to-row 
transfer matrices can be reduced from dimension N L to dimension 2 m , where m is not 
greater than L. Like the Ising model, the partition function is a matrix element of a 
direct product of m matrices, each of dimension two. 

For the Ising model, one can calculate the correlations quite explicity as determi- 
nants by using free-fermion operators. [191 |2"D] or equivalently by writing the partition 
functions directly as pfaffians.|21| 

If the partition functions of the JV-state superintegrable model resemble those of 
the Ising model, then perhaps the correlations are also similar, and can be obtained by 
similar methods. In particular, the order parameter can be defined as 

Mr = Wl, (1-2) 

where r = 1, 2, . . . , TV — 1, 

w = e ' , (1.3) 

and a is a particular spin inside the lattice. We can use this definition for a finite 
lattice: we do so herein. We only expect the simple result (|l.ip to be true in the limit 
when the lattice is large and the spin a is deep inside it. It should then be independent 
of the values of the rapidity parameters, so should be the same for the general solvable 
model as for the superintegrable case. 

We take the spins to have the values 0, . . . , N — 1 and a = {oi, . . . , ox} to be the 
set of spins in a horizontal row of the lattice. Let uj, be the iV L -dimensional vector 
with entries 

(u b )„ = S(a 1 ,b)8(a 2 ,b)---S(a L ,b) (1.4) 
and S r the diagonal matrix with entries 

L 

(S r ) a , a , = ^IpKX). (1.5) 

3=1 

Then with these boundary conditions, (|1.2[) can be written in terms of the row-to-row 
transfer matrix T as 

Mr = 4™^, (1 . 6) 

u\ T m + n ii 

where m (n) is the number of rows below (above) the particular spin a. We have chosen 
a to lie in the first column: since we are using cylindrical (cyclic) boundary conditions, 
this is no restriction. 

The tranfer matrix T commutes with a hamiltonian H. For convenience, we replace 
(fj~6)) with 

M r = ^ ^ ^. (1.7) 

u T e -( a + /3)H Uo 

In the limit of m, n large, the only eigenvectors of T entering the RHS of (|1.6p are 
those corresponding to the iV asymptotically degenerate largest eigenvalues. The same 
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is true of (|1.7p in the limit of a, /3 large and positive. Hence these limits of (|1.6[) , (11.71) 
must be the same. 

Let R be the operator that increases every spin in a row by one: its elements are 

L 
3=1 

Since R N = 1, its eigenvalues are l,u,u 2 , . . . ,u . Let Vp (for P = 0, 1, . . . , TV — 1) 
be the space of vectors v such that 

Rv = lu P v. (1.9) 
Then the full N L -dimensional space is the union of Vo, . . . , Vjv-i- Let vp be the vector 



vp 



N- 1 ' 2 Y J ^ Pb u b . (1.10) 

6=0 

Then vp and e~ l3n vp are in Vp. If 

Q = P + r , mod TV, (1.11) 

then S r e~ /3n vq is also in Vp. Vectors in different spaces Vp, Vq are orthogonal. It 
follows that, for b = 0, ... N - 1, 

JV-1 

^e- aH 5 r e _,3K uo = TV -1 cj _pb Wp,p+ r (a, /3) (1.12a) 

p=o 

JV-1 

N~ 1 ^2lo~ pi ' Zp(a) , (1.12b) 



t -aH 



where 



W PQ (a,l3) = v%B- aH S T %- pH v Q , (1.13a) 



Z P (a) = v r p e~ aH v P . (1.13b) 

Each LHS of (1.12) is the partition function of a lattice where the top spins are fixed 
to the value 0, the bottom to the value b. If b ^ this ensures that there is at least 
one mis-matched seam (between phases where most spins are zero and most spins have 
value b) running horizontally across the lattice. If £ is the interfacial tension, (which 
we expect to be independent of L) then this will make each partition function smaller 
then that for b = by a factor e" ic .[22l section 7.10] 

In the limit of L large the ratio of these expressions for b ^ to their values for 
6 = will therefore become zero. From (1.12), it follows that in this limit Zp(a) is 
independent of P, while Wpq{oi, /3) depends on P, Q only via their difference r — Q — P. 
We show at the end of section 5 that these assertions are certainly true in the limit 
when a, f3,L all tend to infinity. 

The numerator in (|1.7|l can therefore be replaced by Wpq(o, /3), for any P, Q sat- 
isfying HI. lip . The denominator can be replaced by Zp(a + /3), for any P, but if a is 
also large, each Zp(a) is of the form Ke va , where K, v must be independent of P, so 
we can more symmetrically replace the denominator by [Zp(2a)Zg(2l3)] 1 ^ 2 , giving as 
our final expression for the spontaneous magnetization 

M {2) = grgfag) (1 14) 

The three expressions M r , A^' 1 ', M?^ are equal in the limit L,m,n or L,a,j3 all 
becoming infinite. 
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Previous calculations for finite L, a, (3 

In [231 1241 [25] . we have looked at the problem of calculating Wpq{ol, /3), Zp(a), and 
hence Mr , M r 2 \ algebraically, for the superintegrable chiral Potts model, with finite 
L, q, p. We refer to these papers as I, II, III, and prefix their equations accordingly. 

The calculation of Zp(a) is straightforward, being a minor adaptation of the par- 
tition function calculations of [HE], and is given in paper II. The real problem is to 
calculate Wpq(o, /3), or equivalently the ratio 



T>P Q (a,P) = 



Zp( a )Z Q (f3) 



If we also define 



then (|1.14|1 becomes 



Z P (a) = Zp(2a)/Zp(a>y 



Mi 2) = 



T>p Q (a,f3) 



(1.15) 



(1.16) 



(1.17) 



In I we considered the N = 2 case 
spinor operators (Clifford algebra) [26] 



which is the Ising model. We used Kaufman's 
to first write V PQ {a, f3) (for P = 0, Q = r = 1) 
as the square root of an L by L determinant in 1.4.59. Then in section 6 of I, eqn. 

I. 6.29 and 1.7.9, we further reduced this result to an m by m determinant (with no 
square root), where m < L/2. With obvious modifications of notation to allow for the 
working of the later papers, and taking p = 0, x = 1 in 1.3.5, 1.3.6 (also in II. 5. 23 and 

II. 5. 25), this result can be written as 



V PQ (a,f3) = det[/ m - Xp(a)Ep Q BpQXQ((3)EQpB Q p} 



(1.18) 



where I m is the identity matrix of dimension m, Xp(a), Epq, Xq(/3), Eqp are diagonal 



matrices, Bpq is an m by m' matrix, where \m- 
is orthogonal in the sense that 



= or 1, Bqp = -B P q and B P q 



BpqBpq = I„ 



if m > m . 



BpqB P q = /„ 



if m < m' . 



(1.19) 



We used the first (L by L) form to take the limit a, /3 — > +oo: the result of course 
agreed with that of Yang 20 and Montroll, Potts and Ward [21] . and with (|1.1[) above. 

In II we considered the superintegrable chiral Potts model and showed that the 
^^-dimensional matrices in (|1.13a|l could be replaced by ones of lower dimension. In 
particular, the H in the first exponential could be reduced to dimension 2 m , where (for 
P = 0,...,JV-1), 



m — mp = integer part of 



(N - 1)L 



(1.20) 



Similarly, the second H could be replaced by one of dimension 2 m , where m' = niQ, 
and iSy by a 2 m by 2 m matrix Spq. (The p,q,S^ ed of paper II became P,Q,Spq in 
paper III and herein.) 

We went on to conjecture that ()1.18[) . (|1-19|) also applied to the superintegrable 
model, with fairly obvious generalizations of the definitions of the X, E, B matrices. 
We observed that this conjecture agreed with numerical tests performed to 60 digits of 
accuracy. 

These calculations involved sets of m quantities 9\,...,0 m defined by 



cos 9j 



(l + w j )/(l-w j ), 0<6>i<7r, 



(1.21) 
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where Wi, ■ ■ ■ , w m are the zeros of the m-th degree polynomial p P (w) given by 



N-l , N xL 

71 = V 7 

taking w = z . Let c = (1 + w)/(l — tu) and, for all complex numbers c, 

Pp(c) = N- L (c+l) m p p (w). (1.23) 
Then Vp(c) is the polynomial with zeros ci, . . . , c m , i.e. 

m 

Vp{c) = Y[(c-cj). (1.24) 

2=1 

Similarly, we can define m' quantities #i , . . . , 0^, and c[ , . . . , c^, by replacing P by Q 
in the above three equations. The p p (w), Vp(c) here are those of paper III, which are 
the P(w), P p (c) of II. 2. 17, II. 2. 18 and II.6.4. 

In III we showed that both the S r and Spq matrices satisfied various commutation 
relations with the hamiltonians, in particular that Spq satisfied III. 3. 39 and III. 3. 40. 
We conjectured in III. 3. 45 that the elements of Spq were simple ratios of products of 
trigonometric functions of the 6t. We suggested that this result applied for any values 
of the 8i and 6[, not necessarily those given by (|1.21|) and (|1.22[) . Further, if we defined 
Vpq by III. 3. 48, then it was also given as a determinant by III. 4. 9 and III. 4. 10, again 
for arbitrary 6i, 6[. 

Outstanding problems and progress 

It therefore appears that there is indeed an algebraic route to calculating T>pq. How- 
ever, there are still three outstanding problems to be overcome: 

1. To prove that the elements of Spq are given by III. 3. 45, and hence T>pq by 
III.3.48. 

2. To further prove that T>pq is given as a determinant by III. 4. 9 or equivalently 
III.4.10. 

3. To calculate the determinant III. 4. 9 in the limit L, a, f3 —> oo so as to regain the 
known result (|l.ip . This has not previously been done directly even for the N = 2 Ising 
case: in paper I we calculated (|1.1[) from the expression for T>pq as the square root of 
an L by L determinant, using Szego's theorem. This theorem was derived in reponse 
to the first (unpublished) derivation of the Ising model spontaneous magnetization by 
Onsager and Kaufman. 27, 28 It was later used my Montroll, Potts and Ward. [21] 

Progress has been made. We have proved that the expression III. 3. 45 for Spq 
satisfies the commutation relations III. 3. 39 - III. 3. 41. From numerical calulations for 
small N, L (N, L < 6), it appears that these relations (which are linear in the elements 
of Spq , with many more equations than unknowns) determine Spq uniquely. If so, 
then III. 3. 45 and III. 3. 48 have to be correct. 

In paper III we defined yi, y[ to be the elements of the diagonal m and m'-dimensional 
matrices Y = Xp(a)EpQ, Y' = Xq(/3)Eqp. Both the expression III. 3. 48 for Vpq as a 
2 m+m -dimensional sum, and the expression III. 4. 9 as an m-dimensional determinant, 
are rational functions of the a, c' t , yi, y'i- One can take all these variables to be arbi- 
trary and can verify that the denominators are identical. One can then prove that 
the numerators are also the same by a recursive method using the symmetries and the 
fact that if c m = c' m , , then each expression simplifies to one with m, ml replaced by 
m— 1, m'— 1. 

In some ways the hardest of the three problems is to take the limit L,a, ft — >• oo. The 
determinant for Dpq is a hugely smaller calculation than the original 2 L -dimensional 
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sums in (jl . 13a|) . (|1.13b[) . but it is still ultimately infinite. We have succeeded in calcu- 
lating the determinant in the limit a, /3 — > oo as a simple product, for finite L,m,m': 
the key trick is to note that when a, j3 — » oo the matrix sum in (|1.18[) can be written 
as the product of two Cauchy-like matrices. The result follows by then taking the 
L — v oo limit of the product. It is this calculation we report here. We hope to publish 
the work on the first two problems later. 



2 The matrices 

We shall need the definitions of the X and E diagonal matrices. From II. 3. 16 and II. 7. 4 

\v I W _ -fc'sinflj sinh(iVaAi) 8 uj 

1 Pl J|,J A l cosh(iVaA i ) + (l-fc'cos6» I )sinh(iVQA l ) ' ( ' 

where 

\i = (l-2fc'cos6» 1 + fe' 2 ) 1/2 . (2.2) 
The matrix Xq(P) is defined similarly, with P, a, 6i, \i replaced by Q, /3, 6' iy Xj. 

From II. 6. 18 and II. 6. 19, the matrix Epq is an m by m diagonal matrix with entries 

[Epgh = e(P y Q,i)S 1 , J , (2.3) 
where, for < P, Q < N and P / Q, 

e(P,Q,i)= sin Si if P < Q and m' = m — 1 

= tan((9i/2) if P < Q and m' = m 

= l/sin0i if P > Q and m = m+1 (2.4) 

= cot(0i/2) if P > Q and m! = m. 

These equations cover all cases. The matrix E(Q, P) is defined similarly, with P, Q 
interchanged, m,m! also interchanged, and 9i replaced by 9'i. 
Bpq is an m by m Cauchy-like matrix with elements 

(Bp Q h = -^K- (2.5) 

a- Cj 

Given ci, . . . , c m , c'i, . . . , c' m , with \m — m'\ < 1, there is a unique way of choosing 
/i, . . . , fm, f[, ■ ■ ■ , f' m i so that Bpq satisfies the orthogonality condition (|1.19[) . The 
working is given in section 6 of II. We remark in section 4 of III that it is true for 
arbitrary Ci , c[ . Let 

m m 

a« = JJ(Cf - Cj-) , a't = J^[( c » - c j) ) ( 2 -6) 

3=1 3=1 
m m' 

h= Y[ ( c »~ c j)' b 'i = II ( c »- c j)' 

then the results II. 6. 8, II. 6. 13, II. 6. 16 can be written as 

fi = ea t /bi , f- 2 = -ea'i/b'i, (2.7) 

where e = ±1 is independent of i. 

( For the particular values of d , c[ given by (|1.2ip and (|1.22[1 , we observe numerically 
that and / 8 ' 2 are positive real if we choose e = 1 if P < Q, and e = — lifP>Q.) 

A quantity that we shall need is 

A f„ >\ _ ni<i<j<m( C ' " C j)IIl<i<j<m'( C 3- ~ gf) , v 
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3 The function Zp(a) 

The partition function Zp(pt) is, from II. 3. 16 and II. 5.38, or from III. 3. 27 and III. 3. 29, 

ry , \ -» P oTT Ai coshfiVaAi) + (1 — fe' cos sinh(iVaAi) ,„ lN 

Z P {a) = e j| , (3.1) 

i=l ' 

where 

/i P = 2fc'P+(l + fc')(miV-iVL + L). (3.2) 

When a is large, Zp(a) has the form Ce ga , where C, g are independent of a and 

m 

S = gp = -/ip + X! Al ' ( 3 ' 3 ) 

i=l 

Hence from (I1.16|) . Zp(a) — > C~ 1 as a — > oo. 
Set 

Xp = Xp(oo), ii = (lp)i,i , 2p = 2p(oo) = C- 1 (3.4) 
and let c, A be variables related to one another by 

c = cos 6» = (1 + fc' 2 - A 2 )/2fc', (3.5) 

then the a, 0i, Xi of (11.211) . (12. 2[) are related in the same manner. Instead of viewing Xi, 
-Zp, etc. as two- valued functions of the Cj, we can regard them as single- valued rational 
functions of the Ai. Then 

2 _ fc' 2 -(l-A0 2 



(l + A0 2 -fc' 2 ' 
4 A, 



(3.6) 



^ = n (1+ ^. fc , 3 = n<i+*?>- 

i=l v ' i=l 

Analogous relations apply, with P,\i,Xi replaced by Q,Ai,a4, respectively. 
Another function that we shall useful is 

together with the elementary identity 



^m,m'(A 2 ,A' ) 



,(A,A') 



> 



-.(m-m') 2 TT 



n^ 2 rn^m) . <w 

»=i * / j=i 



4 Calculation of T>pq 

For any m by m' matrix A, and m' by m matrix B, it is true that 

det(J m + AB) = det(J m / + BA) , (4.1) 

so from (fl~T8)l . 

Dp Q (a,^) = det[J m /-Jfg(/3)SQp5 Q pJfp(o)SpgBpQ] 

= VqpOS.o). (4.2) 

This symmetry also follows directly from the definitions (|1.5|) - (|1.15[) . the fact that H 
is hermitian and St — 5_ r . 
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Without loss of generality, we can therefore restrict our attention to the case P > Q, 
when m < m' . 

Then from (|1.19[1 we can write I m in (|1.18[) as BpqBJ,q. Remembering that Bqp — 
—BpQ, we can then write (|1.18[) as 

V PQ (a,l3) = det[UBl Q ], (4.3) 

where 

U = B PQ + Xp(a)Ep Q Bp Q X Q (P)E Q p . (4.4) 

Define 

y t = [X P (a)]^e(P,Q,i), y< = [X Q (p)]jj e(Q,P, j) , (4.5) 
then from (12.51) the elements of U are 



Utj = f*m+V^ m ( 4. 6 ) 

Ci Cj 

In general we do not know how to calculate the determinant of such a matrix. 
However, if we take the limits a, /3 — > +oo and express Ci,c'j,yi,y'j as rational functions 
of Xi, we find that a factor A 4 + X'j cancels out of the RHS of (|4.6[) . If m' — m, the 
result is Cauchy-like matrix, and one can calculate the determinant of U. 

Hereinafter we take the limit a, /3 — i- +oo, so 



Xie(P,Q,i), y'j = x'j e(Q, P, j) , (4.7) 



where Xi, x'j are given by (|3.4p . We write T>pq(x,, oo) simply as T>pq. The integer L is 
still finite. 

The case P > Q, m = w! 

The simplest case is when m — m' and all matrices are square, so from (|1.19|1 and 
(E3J), 

Vpq = detU/ det Bpg . (4.8) 

Cauchy-like matrices 

If A is the mbyra matrix with entries 

A %] = — l —, (4.9) 

then it is a Cauchy matrix and its determinant is A m , m (c, c'), using the definition 
(|2.8p .[2l?l eq. 2.7] Any matrix with elements of the form ()2.5|l is said to be Cauchy-like, 
and has determinant 

m 

detBpg = A m , m (c,c) Y[fifi (4.10) 

i=i 

for all fi,fi- We have in fact chosen the fi,fi so that Bpq is orthogonal, so has 
determinant ±1. However, the form (|4.10p is convenient here as the /j, products will 
cancel out of (|4.8[) . 
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The determinant of U 



From (|4.8[1 . we still have to calculate the determinant of U. Its elements are given by 
(14. 6[) . so U is not in general Cauchy-like. However, in the limit a, /3 — > oo we find that 
a common factor cancels from the numerator and denominator of (14. 6[) , and U becomes 
Cauchy-like. We can then evaluate its determinant by parallelling (|4.9|l - (|4.10l) . 
The case P > Q, m — m' is the fourth one listed in (|2,4[l . so 

(Ep Q )ij = cot(6»i/2) <5y , (S QP )y = tan(0-/2) ^ . (4.11) 

From (|2.ip . taking the limit a — > +oo, 

-l«=rM^' (4 ' 12) 

Noting that sinfli cot(0»/2) = l + a = [(l + k') 2 - Xi 2 ]/2k' , it follows from (gjj) that 

1 + k' - A, 

w = - i-k' + Xt - (4 ' 13) 

(A factor 1 + k' + Xi has cancelled.) 

The calculation of y'j is similar, except that now we use sin^ tan(#J/2) = 1 — c'j = 
[Xj 2 - (1 - k') 2 ]/2k' to obtain 

* = rrvri • (4 - 14) 

Thus 

1 + ^ = (i-.^+Ua;) - (4 - 15) 

Also, 

= (Aj- 2 — X 2 )/2k' . (4.16) 
We see that the factor Ai + Aj cancels out of (|4.6|l , leaving 

u.. = ~ 4k 'f^ um 

3 (l-fc' + A l )(l + fe' + A;)(A l -A;) { ' 

so {/ is a Cauchy-like matrix, similar to Bpq, but with the denominator Ci — c'j replaced 
by A; — X'j. Analogously to (|4.1Up . its determinant is 

dett/ = A_(A,A') 11 (1-^ + ^(1^ + At) ' (4J8) 

Z= 1 



Also, from (|4.16p , we can write (|2.5|) as 

A? -A 



(Bp Q h = (4.19) 



j 

so det Bpq is also equal to 



det Bpq = A m , m (A 2 ,A' 2 ) jj(-2fc' /</i) • (4-20) 



The /i, /j cancel out of the ratio (|4.8I) . leaving 



T>pq = 



(A, A') 



#\ II (i _ k' + Xi)(l + k' + A') ■ (4 - 2i) 



, m (A 2 , A' ) 11 (1 - fe' + Ai)(l + k' + AO 

i= 1 

From (fi~17| . ((3T7ll . (jH| and (H) it follows that 

(tm (2) ) 2 - ^±^n£i^j (4 22) 
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All cases 



When P > Q and m' = m+1, we can still write 1>pq as in (|4.3[) . However, [/ and 
Bpq are no longer square matrices, so we can longer simply take products or ratios of 
determinants, as in ()4.8|) . Even so, we can still calculate T>pq (and hence M^) by 
adding a row to Bpq and U to make them square Cauchy-like matrices. The matrix 
UBpQ in (|4.3[) is then m' by m', but is of upper- block-triangular form. The top left 
block is the original mbjra matrix UBpQ, while the lower-right block is the one by 
one unit matrix. Hence the determinant (|4.3p is unchanged and can be evaluated as a 
product of the two m' by ml determinants. 

We do this in the Appendix. The result (|A18|) is the same as (|4.22[) . except that 
the factor 72.(1 + k') is inverted. 

From (|1.17|) and (|4.2[) , is unchanged by interchanging P with Q, m with m' 

and the \i with the Xj. From (|3.8[) . this inverts the function 1Z. We can use this 
symmetry to calculate M.^ in the other two cases. 

For the four cases (|2.4|) . define the factor Q by 



g = n(i-k')n(i + k') 


if 


P < Q 


and 


m 


= m— 1 


= 11(1 - k')/K(l + k') 


if 


P <Q 


and 


in 


= m 


= l/[fc(l - k')1l(l + k')] 


if 


P>Q 


and 


/ 

m 


= m+1 


= 11(1 + k')/Tl{l - k') 


if 


P>Q 


and 


i 

m 


= m , 



then we find that 

m / m 

(m^) 2 = gjinx'j) Xlnx) (4.24) 

3=1 I i=l 

for all four cases. 



5 The limit L — >• oo 

The result (|4.24|l is exact for finite L. The last step in the calculation is to let L — > oo. 
For this we shall need for the first time herein the particular definition (|1.22[) - (|1.23[) 
of the polynomial (|1.24|l . 

Let c, A be two variables related as are Ci,Xi in (|3.5[1 . Noting that — A 2 = 
2k' (c - a), it follows from ((3~8|l and (fl~24)l that 

H(\)K(-\) = (k'/2r' m 'Pp(c)/P Q (c). (5.1) 

Using (|1.22p . p.23p . this can be written as 

Tl(\)TZ(-\) = [fc'(c+l)/2] m - m ' z Q ' p yV{\), (5.2) 

where 

W(A) = Wp{X)/Wq(\), (5.3) 
jv-i L 

wp(A) = i+E- nCF+i) (^) ■ ( 5 - 4 ) 

n=l 

We have cancelled the n = term in (|1.22[) from the ratio (15. 3[) . 

These z, w, c are related to A by the relations (|1.22[) . (|1.23p . (I3.5[) . In particular, 

z - w - c + i - (i + k>y-x>- (5 ' 5) 
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When A 2 < (1 — k') 2 we choose z to be positive real. Thus w,c are rational functions 
of A. We see that z is multi-valued, but we can choose it to be analytic by cutting the 
complex A-plane from 1 — k' to 1 + k', and from — 1 — k' to — 1 + k' ■ It is then analytic 
on the imaginary axis and at infinity. 

The function W(A) is even and rational, its poles and zeros being symmetrical about 
the imaginary axis, with none on the axis. Also, z — > 1 and W(A) -> 1 as A -> oo, so 
log W(A) is analytic in a vertical strip containing the imaginary axis, and tends to zero 
at infinity as 1/A 2 . 

We can therefore perform a Wiener-Hopf factorization[30] of W(A): 



W(A) = W+(A)W_(A) : 



where, for Sft(A) > 0, 



logW+(A) = - 



2tti 



logW(A) JA/ 
A' -A dA 



(5.6) 



(5.7) 



the integration being up a vertical line just to the left of the imaginary axis. This 
W+ (A) is analytic and non-zero in the right-half complex A-plane and on the imaginary 
axis, and tends to one as A — > oo. The function W_(A) is defined similarly, but with 
3?(A) < and the sign of the RHS reversed: it is analytic and non-zero in the LHP and 
on the imaginary axis. 

From 1)3. 8[) . the function 72(A) is analytic in the RHP and is proportional to (A/2) m ~ m 
when A — > oo. It follows from (|5.5[) and 

(1 + fc') 2 — A 2 



that 



72(A) 



fc'(c + l)/2 = 
1+A+fc' 



1 - k' + A 
1 + k 1 + A 



(Q-PJ/JV 



W+(A). 



(5.8) 



This is an exact result, true for finite L. 
Let 



t 



k' 



2/JV 



1 + k' 

then when A lies on the imaginary axis, z lies on the positive real axis, between t and 
1. Hence, for < n < N, 

i - ■ i _ /- 

<1. (5.9) 



1 - z 


< 


1 


- 1 


LU n — Z 


0J 


- 1 



It follows from (|5.3[l . (15. 4[) that W(A) tends uniformly to one as L — > oo. For sufficiently 
small e, this must also be true for A on the integration line in (|5.7[) . Hence 



W+(A) ->■ 1 as L ->■ oo , for Sf(A) > , 



so ^(A) is then given by (JS3J with W+(A) = 1. 
It follows that 



72.(1 + &' 



/\1 {Q-p)/n 



TZ(1 



72.(1 + fe') = ( 1 + fc ' ) m-»'+(F-Q)/JV ) ^!_ fc /) = (i_fc') (Q " P)/JV . 
Using r = Q — P, modulo JV, we obtain 

fc /2-.r(JV-r)/iV 2 



(M 2) ) 2 = (1 



(5.10) 

(5.11) 
(5.12) 

(5.13) 



for all four cases, in agreement with (11. ip . 
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Asymptotic degeneracy 

We are now in a position to confirm the remarks we made before (|1.14[) . Expanding 
(|3.8[) and (|5.8|) to first order in 1/A for A large and equating the coefficients, using (|3.2[) 
and (13.311. we obtain 



gp = 9Q ■ 

We noted after (|3.3|) that when a is large, then Zp(a) is proportional to exp(gpa), 
so Zp(a) I Zq(ol) must tend to a limit as a — > oo. From (|1.16|l and (|3.4|) . this is Zq/Zp, 
so from (|3~71) and (f378l) , 



w = «ft-W + »0 = (514) 

using (|5,8[l with W+(A) = 1 to evaluate the middle expression. 

Hence in the limit a,f3,L — >• oo, Zp(ot) is the same for all P. We have already 
shown in (|5.13[) that depends on P and Q only via the difference r. From (|1.15[) 

- ()1.17[) . the same must be true of Wpq(oi, 0). This justifies our replacing the sums in 
(1.12) (for b = 0) by single terms. 



6 Summary 

The magnetization Air of the chiral Potts model can be expressed as the expectation 
value of uj ra , where a is a spin inside a cylindrical lattice with the spins at the top 
and bottom boundaries fixed to zero. This was calculated in [3l S] analytically by 
generalizing M r , showing that it satisfied functional relations in the large-lattice limit, 
and then solving those relations. This was quite a different method from the algebraic 
techniques originally used for the Ising model. 

To obtain M r for the general solvable model it is sufficient to obtain it for the 
superintegrable case, and this case has many resemblances to the Ising model. This 
naturally leads to the question whether there is an algebraic way of calculating the 
magnetization for the superintegrable chiral Potts model. 

We looked at this problem in three previous papers, 23 - 25] which we refer to as I, 

II, III. In I we revisited the Ising model (the N — 2 case of the superintegrable model) 
and showed that M r was proportional to the determinant 1.7.7 and 1.7.9. 

In II we first showed in II. 5. 37 that M r was proportional to the weighted sum of 
the elements of a 2 m by 2 m matrix S^ cd . Generalizing our calculation in I, we then 
conjectured that it was proportional to the m by m' determinant T>pq(ol, ft) of (|1.18p 
herein. 

In III we further showed that the matrix Spq = S^ ed satisfied a number of com- 
mutation relations, in particular the two relations (2.21), (2.22) therein, and further 
conjectured that the elements of Spq had the particular simple product form III. 3. 45. 

We have since proved that this form for Spq both satisfies the two commutation 
relations and implies the determinantal result. We hope to publish the working soon. 
What we have not done, and would be needed to complete the proof, is to show that 
the commutation relations III. 3. 39 and III. 3. 40 (plus the simple normalization property 

III. 3. 41) define Spq uniquely. Calculations for small lattices suggest that this is so. 
These calculations are all for finite sums and determinants of finite matrices. They 

involve parameters a, f3, which can be regarded as a measure of the number of rows 
below and above the selected spin a in the lattice. They also involve the lattice width 
L. Ultimately we want to take the limits a, /3, L — > oo. 

To complete this algebraic calculation of Ai r , we need to evaluate the determinant 
T>pq(o., ft). We do not know how to do this for finite a, ft, but we show in section 4 
and the Appendix herein that in the limit when a and /3 are both infinite, T>pq(q., j3) 
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can be written as the determinant of the product of two square Cauchy-like matrices. 
We can therefore evaluate the determinant as a product of terms, the number of terms 
being quadratic in m,m! . Finally in section 5 we take the limit L — > oo, showing that 
the needed function 7Z(X) then has a simple form for 5ft(A) > 0. We of course regain 
the result of [3] [4], which had been conjectured by Albertini et al in 1989. [T] 

For the N = 2 Ising case, this appears to be a new way of calculating the needed de- 
terminant. Yang 20 did so by calculating the eigenvalues in the large- L limit. Montroll 
et aZ,[2T] and presumably Onsager and Kaufman in 1949. [27l 128] did so by expressing 
the result in terms of a Toeplitz determinant and then using Szego's theorem. 131) 



7 Additions 

Since posting this paper on the Los Alamos archive, there has been considerable further 
progress on the first two problems listed above in the Introduction. Iorgov et al 32 have 
proved that Spq is indeed given by III. 4. 9, and therefore Dpq by III. 3. 48. They did 
this by showing that III. 3. 45 satisfied the commutation relations III. 3. 39 - III. 3. 41, and 
were able to prove that these relations have a unique solution. Further, they went on 
to calculate Dpq in the limit a, /3 — > oo directly, using the expression III. 3. 48 as a sum 
over matrix elements, rather than the determinantal form III. 4. 9 or III. 4. 10. 

The author has also posted a paper on the archive 33 giving the proof of the equiv- 
alence of the sum and determinantal forms of Dpq that is outlined in the penultimate 
paragraph of the introduction above. It is remarked therein that the motivation for 
this work was the similarity between the Ising and superintegrable chiral Potts models, 
in particular to find a derivation of the determinant form of Dpq. This determinantal 
form was conjectured in II by generalizing the Ising model result as expressed in I. 
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Appendix A 

Here we consider the case when P > Q and m' = m + 1, so the matrices Bpq, U are 
not square. Then from (|2.4j) . e(P,Q,i) = l/sin^ and e(Q,P,j) = sin#j. Using (I4.5|) . 
(14.121) . we obtain 



-2k' , _ (1 — \'j) 2 — k' 2 

(1 + \i) 2 - k' 2 ' Vj ~ 2k' 



= . , . , ::; , - »', = " '1, ■ m 

Hence 



1 + (l + A,)'-*" ■ (A2) 
Using (|4.16[) . we see that the factor Xi + X'j again cancels in (|4.6[) . leaving 

U = U w + U {2) , (A3) 
where U W ,U {2) are matrices with elements 

rjW = - 4fc 7*/j fA4) 

lJ (A i -A')[(l + A l ) 2 -fc' 2 )] ' V ; 
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(l + A0 2 -fc' 2 ) 
respectively. Thus 

U {2) = £r/ T , (A6) 

where £,rj are vectors and (rj)j — /j. 

Consider the vector Bpq t\: from (|2.5[) it has entries 



(Bp Q v)i = fiT'(ci), 

where 

™ - T,^ 

3=1 3 

Remembering that P > Q, this is precisely the sum considered in (II. 6. 6), but with p, q 
therein replaced by Q,P, so (also using II. 6. 14), 

JF'(c) = 1 ' + V P (c)/V Q (c), (A7) 

where 7' is independent of c. Taking the limit c — > 00, we obtain 7' = 0. It follows 
that J-'(c) vanishes when c = Ci (i = 1, . . . , m), so 

Bpqv = 0. (A8) 

Substituting ()A3|) into (|4.3[1 . the f/' 2 ' term is zero, leaving 

V PQ = det[U w B T PQ \. (A9) 

We have reduced the problem to one of calculating the determinant of a product of 
two Cauchy-like matrices, but unfortunately they are not square. The solution to this 
problem is actually suggested by ()A8|I : of we define an m! by m! matrix 

s - ( y ) , 

(dropping the suffixes P, Q). Then, using 1)1.190 . all its m! rows are mutually orthogonal. 
Multiplying (|A7)I by c and taking the limit c — » 00, we obtain 

m 

j=i 

so together with ()A8[I it follows that BB T = I m >, i.e. B is a square orthogonal matrix. 
We also extend U^ 1 ' by adding the row r/ T to form the square matrix 

"-[*)■ (A12) 

Then, using lfM| and (1A11)) . 

UB r _ ( U"Bl a BO), \ _ / BWIS, C™,, \ 

V >7 T S?0 » T 1 / \ 1 / 

We see that UB T is an upper-right block triangular matrix and det U B T = det U^ 1 ' BpQ. 
From (|A9P , using the orthogonality of B, 

V PQ = det Uj det B. (A14) 
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The square matrices B and U are Cauchy-like. All the elements Bij of B are given 
by the RHS of (|4.19|) . provided we take f m +i = — X 2 n+1 /2k' and then let A m +i — > oo. 
Using the general formula (|4.10[) for an m! by m! determinant, and then taking this 
limit, we obtain 

m m 

det B = A m , m , (A 2 , A' 2 ) [[(2*70 JJ f s , (A15) 

i=l 3=1 

where A m>m /(c, c') is defined by (|2,8|) ■ (Similarly to section4, we know from the or- 
thogonality of B that its determinant in ±1, but this form is convenient here because 
the fi, fj products will cancel out of (|A14[) .~) 

Similarly, all the elements Uij of hi are given by the RHS of (|A4[) . except that now 
we take f m +i = — A^, +1 /4fc' before letting A m +i — > oo. Again using the general formula 
P~T0)l . we find that 



-k' 2 

3=1 

Hence from (|A14I) . 



detW = A m , m ,(A,A') J] (TT0^ Iltf- 



(A16) 



27 



PQ 



<(A,A') 



iA--.A-o n (1 + Ai)2 „ fc ,2- (Air) 



Using (TrT7)) . (J37J, (|3^8)l and p^|l, we now obtain 



(A^n 2 nr=i^( A i) rA18 x 

1 r) n(i +k>) -k>) ^^(Ai) ■ 1 J 



This is the result quoted at the end of section 4. 
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